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^ • Abstract 

^ • We report a strain-controlled tuning of magnetism in transition-metal-atom-decorated graphene. 
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Our first-principles calculations demonstrate that strain can lead to a sudden change in the mag- 
netic configuration of a transition metal (TM) adatom and the local atomic structure in the sur- 



'^ ■ rounding graphene layer, which have a dramatic effect on the effective exchange coupling between 

O . neighboring TM atoms. A strong spin-dependent hybridization between TM d and graphene vr 



orbital states, derived from the orbital selection rule of the local lattice symmetry, is responsible 



^ ' for the determination of the local electronic and magnetic structure. Our results indicate that the 

^^ • strain can be an effective way to control the magnetism of atomic-scale nanostructures, where the 

'^ ! reliable control of their magnetic states is a key step for the future spintronic applications. 
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A reliable control of magnetic states is central to the use of magnetic nanostructures in 
future spintronics and quantum information devices IL Although it is shown that magnetism 
could be generally modulated by external magnetic iHs] or electric j^ fields, it is desirable 
to find an alternative scheme to control the magnetism for various spintronic applications. 
Recently, graphene has gathered tremendous attention due to its unique electronic and 
mechanical properties for nanoscale electronics [s!]. As a candidate material for spintronic 
devices, transition-metal-atom-decorated graphene (denoted as TM-graphene hereafter) has 






been studied extensively in theory p-lsl and experiment |9l-lll| . manifesting some remarkable 
electronic and magnetic behaviors. Thus, developing a novel method to tune the magnetism 
of TM-graphene system is quite urgent for future spintronics applications. 

Graphene is the thinnest material ever synthesized [Sj]. While it is one of the strongest 
materials ever measured experimentally 12|-|l5|. it has been shown that it can sustain elastic 
deformations as large as 25%. In this Letter, we predict that the strain is an effective way 
to control the magnetic properties of TM-graphene systems. As for the graphene layers, we 
considered a pristine graphene (PG) layer as well as graphene layers with defects of a single 
vacancy (SV) or a double vacancy (DV)- It has been shown that vacancies in graphene can 
be created by electron irradiation 16|, [iTJ] or ion bombardment 18j. Our results for the Mn- 
atom-decorated graphene demonstrate that the strain can induce a sudden change of the 
local atomic structure of graphene around a Mn atom and the spin state of Mn atom, which 
can lead to a striking change of the effective exchange coupling between neighboring Mn 
magnetic moments. A strong spin-dependent hybridization between Mn d and graphene tt 
orbital or dangling bond states is responsible for the determination of the local electronic and 
magnetic structure. We also showed that the strain-controlled magnetism in Mn-graphene 
is a general phenomenon happening in different TM-graphene systems. 

All the density-functional-theory (DFT) calculations were performed using VASP 
code'l9]. We used the projector augmented wave (PAW) potentials and the generalized 
gradient approximation with the Perdew-Burk-Ernzerhof (PBE) functional to describe the 



core electrons and the exchange-corre 



ation energy, respectively, which were proved to de- 



scribe well the TM-graphene systems 7|, |20|, |2l|]. A 7x7 graphene supercell (98 C atoms) 
was used in our calculations, approaching the single defect hmit ( 1% defect concentra- 
tion). A 12x12 supercell was used to calculate the interaction of TM atoms on graphene. 
A P-centered 6x6x1 k-point sampling was used for Brillouin-zone integration. Moreover 
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FIG. 1: (a) Magnetic Moments of Mn atoms adsorption on prefect graphene (PG) and defected 
graphene with a single-vacancy (SV) as a function of strain, (b) The heights between Mn atoms 
and the nearest C atoms as a function of strain for Mn-PG and Mn-SV systems, (c) The average 
bond lengths between Mn atoms and the nearest C atoms as a function of strain for Mn-PG and 
Mn-SV systems, (d) The binding energies of Mn atoms on graphene as a function of strain for 
Mn-PG and Mn-SV systems, (e) The local structure transition for Mn-SV system from rj = (up 
structure) to rj = -0.04 (down structure). 

GGA+U (with the on-site Hubbard U correction) calculations with U = 5 eV were per- 
formed to confirm the GGA results of the Mn-graphene system. We verified that both GGA 
and GGA+f/ results are consistent with each other. 

We take the case of Mn atom adsorption on graphene as a primary example to demon- 
strate the strain effect on the TM-graphene system. The calculated strain-dependent mag- 
netic moments are shown in Fig. la. A biaxial strain rj is defined as r^ = Aa/a, where a is the 
lattice constant of free standing graphene. A positive (negative) t] means the TM-graphene 
system is under tension (compression). For PG without vacancies, Mn atom prefers to stay 
at the hollow site (above the carbon hexagon center) of PG and the total magnetic moment 
of Mn atom is 5.7 fiB, indicating that Mn atom is in a high-spin state. The height and 



average bond length between Mn atom and its nearest C atoms is 2.09 A and 2.54 A, re- 
spectively, as shown in Fig. lb and Ic. The binding energy of Mn on PG is about 0.45 eV, 
as shown in Fig. Id. When a biaxial tensile strain t] is applied for Mn-PG, interestingly, 
the magnetic ground state changes from high-spin to low-spin (~ 3.1 fis) abruptly at the 
critical strain oi r] = 0.05, as shown in Fig. la. Accompanying the magnetic transition, the 
height and average bond length between Mn atom and its nearest C atoms decrease sharply 
to 1.50 and 2.16 A, respectively, as shown in Fig. lb and Ic. The binding energy of Mn 
atom of Mn-PG increases to 1.15 eV when t] = 0.05. In contrast to the sharp transitions of 
spin-states and local adsorption structures, the binding energy of Mn atom on PG increases 
continuously as a function of rj, as shown in Fig. Id. 

TM atoms on PG have low migration barriers so that it may be mobile at room tern 
perature [7|, [ll|. 
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21| . However, our climbing-image nudged-elastic-band calculations 22| 
show that the migration barrier for the Mn atom on PG increases from 0.36 eV to 0.89 eV 
under a tension oi rj = 0.10, giving rise to a decrease of more than 8 orders of magnitude 
in the diffusion coefficient at room temperature, which implies that the migration can be 
suppressed significantly. It also indicates that an adequate variation of strain can be used to 
control the patterning of TM atoms on graphene. Furthermore, it is known that the mobility 
of the TM atom becomes markedly reduced in the presence of vacancies. For instance, the 
migration barriers for TM atoms on vacancies are larger than 3 eV[7|. 

A similar magnetic transition between high-spin and low-spin states is also found for Mn 
atoms adsorbed near the SV sites in graphene (Mn-SV). Different from the Mn-PG case, 
the unstrained Mn-SV is in a low-spin state with the magnetic moment of ~ 3 /i^, as shown 
in Fig. la. The local adsorption structures of the Mn atoms are shown in Fig. le. The 
binding energy between Mn atom and vacancy is 6.4 eV, which is much larger compared to 
the Mn-PG case, indicating that the Mn-SV system is quite stable. This is in agreement 
with previous results 7|. Interestingly, a compression of r/ = —0.04 can convert the magnetic 
ground state from low-spin (~ 3 /i^) to high-spin (~ 5 /i^) for Mn-SV. Under the compression 
the graphene layer produces a spontaneous rippling with the formation of a C-C dimer bond 
near SV, as shown in Fig. le. As results, the height and the average bond length between 
Mn and its nearest C atoms increase significantly under compression, which weakens the 
binding between Mn atom and vacancy. In general, we find that tensile (compressive) strain 
can enhance (reduce) the binding strength between TM atom and graphene (Fig. Id), which 



C7t f / 



Cnij 



Mnsl 



MniH 



C 7C T ,,+ 



Mnit 



' Mndl 



Cnl 



d'ts't 



d"is"i 



4^ b 



Cti t.. 



Cnlf 



Must 



'%\Mn si 
/ Mr\di 



Ciit/ \ z' yfyindt Cnl 



d't 



dH 



Cti t 




DA„t 



Cji J. 



CJt t 



<e 



Dar^ 



CTI i 



Dc t ,. 



Dei, 



.. Mn di 



W' ^ Mndl 

., z' ny X y' ■ 



DnT 



Dvi, 



y;^j,\\Miirft 

CTtt / \ yz zx ,,:?^^ CTTJ. 



/^=--,\Mnrft ^ , 

' z;i_xv_!r-v '::■ ; L 71 4- 



DA 



c n t V 



Det 



D„i 



dl 



dU 



d'T 



d"i 



FIG. 2: Schematic drawings of the energy diagrams of Mn atom and graphene electronic states 
near the Fermi level (Ef) for (a) Mn atom on PG under rj = 0, (b) Mn atom on PG under rj = 
0.05, (c) Mn atom on graphene with SV under ry = 0, and (d) Mn atom on graphene with SV 
under rj = -0.04. The local electronic structure of graphene is represented by the occupied tt and 
unoccupied vr* levels arising from C p^^ bands, which touch at the K point in the orignal Brillouin 
zone of graphene. In (c), three dangling bond orbitals of C atoms near SV form one bonding state, 
labelled as Da, and two degenerate non-bonding states, labelled as D^. While the bonding Da 
level is close to vr, the non-bonding D^ level is located just above Ep. In (d), the distortion and 
formation of C-C pair gives rise to the bonding and anti-bonding states Db and Dab, which stay 
far away from the Fermi level, and one non-bonding dangling bond state D^ just above Ep. 



is consistent with an experimental observation 11|. 

It is interesting to observe that the electronic and magnetic structures of the Mn atom 
depend strongly on the strain of graphene layer, thereby leading to the strain-induced sudden 
change of the magnetic ground state. The strain affects the local geometry of carbon atoms 
around the Mn atom in two ways. It changes the height of Mn atom above the graphene 
layer and generates a ripple structure under the compressive strain. The change of the local 



environment, which determines the distance and symmetry of C p orbitals with respect to 
that of Mn d orbitals, can lead to a dramatic change in the magnetic d orbitals due to the 
strong spin-dependent hybridization between Mn d and C tt orbital states. 

The isolated Mn atom has s^'^d^^s^^d^^ configuration with a spin exchange splitting around 
4 eV. Since the Mn atom sits over the hollow site, the energy levels of 3(i orbitals are 
broadened by the ligand field from the hybridization between Mn "id and C p^r orbitals of 
the hexagon carbon atoms at the hollow site. Moreover, because the Mn spin-down energy 
level is above the graphene Fermi energy, one electron is transferred from Mn to the graphene 
layer. This is consistent with the observation that the magnetic configuration of the Mn atom 
in the Mn-PG case can be ascribed to d^"^ s^^ d^^ s^^ with the calculated magnetic moment of 
5.7 ;Ub shown in Fig. la. 

In the Mn-PG system, the local geometry of the Mn adatom above the center of the 
carbon hexagon has the Cg^ symmetry. While Mn dyz^zx and d^y^x^-y^ orbitals belong to 
the El and E2 representations of Cq^, respectively, the p^ electrons near the Fermi level 
(Ep) can be also classified to Ai, Ei (corresponding to the vr state) and B2, E2 {tt* state). 
From the orbital selection rule, for example, the C vr* state couples most strongly with Mn 
dxy,x^-y^ orbitals because both of them belong to the same symmetry representation. Due 
to the exchange splitting of Mn d states, however, the spin-up Mn d^ levels are close to the 
occupied C tt states and the spin-down Mn d^ levels are located near the unoccupied C vr* 
states. Consequently, the proximity between Mn c?^ and C tt* states makes the di^ ^2_, 2-71"* 
hybridization much stronger so that the unoccupied C„ ^2_„2 state is pushed down close to 
Ep, as illustrated in Fig. 2a. 

When the graphene is under a tensile strain, the height of Mn atom lowers and the bond 
lengths between Mn and C atoms decreases (Figs, lb and Ic). At r^ ~ 0.05 the change 
becomes abrupt. This is because the reduced bond length with the lower height of Mn 
increases the hybridization strength between Mn d^y^x^-y^^z^ and C vr* states so that the 
dxy,x^-y^ spin-down state becomes occupied. Since the Mn dxy^x^-y^ can hold two electrons, 
consequently, the occupation of the spin-down Mn d^y^x^-y^ orbital reduces the total magnetic 
moment and the exchange splitting as well. Therefore the final magnetic configuration of 
Mn-PG for 1] > 0.05 becomes d^^d"^^ with the magnetic moment of 3 /ib. 

In the case of Mn-SV, one has to consider the role of C dangling bond states arising 
from the carbon vacancy as illustrated in Fig. 2c and 2d. At rj = 0.0, the new structure 



can form without a significant lattice distortion. When a carbon vacancy is present within 
the graphene layer, the Mn atom is much closer to the graphene plane and the C-Mn bond 
length is much shorter. In this case, the main hybridization is between the Mn d state 
and the C dangling bond states from the neighboring C atoms. Similar to the case of 
Mn atom on PG, the coupling between the dangling bond state Da and Df, and the Mn d 
states depends on the symmetries of the states. For example, the largest coupling occurs 
between Mn dyz,zx and C D^ states, whereas Mn d^i couples mostly to Da- Here we have 
a similar spin-dependent hybridization effect due to the proximity of Mn d^ to Da and Mn 
d^ to De levels, respectively. Due to the strong dyz^zx-De coupling, the bonding-antibonding 
separation of the dyzzx-De hybrid states becomes much stronger than in the Mn-PG case, 
so that the antibonding Mn d^^.^^- derived state becomes unoccupied, whereas the bonding 
Mn (i^2 ^x'derived state is occupied as shown in Fig. 2c. Although the occupied bonding 
states are doubly degenerate, its covalent character reduces the effective spin moment close 
to 1 yUB- Thus, as illustrated in Fig. 2c, the magnetic configuration for the case of Mn-SV 
without strain becomes d^^ and d^^ (the Mn s orbital are in higher energy due to the charge 
transfer form Mn to graphene), which corresponds to the calculated magnetic moment of 3 
/iB shown in Fig. la. 

To achieve the high-spin state in Mn-SV, as we have learned from the Mn-PG case, we 
should increase the Mn-C distance to reduce the hybridization between the Mn d and C 
Ptt state. Ironically, this should be done by compressive strain. Indeed, we find that for 
7] < —0.04, the hight of Mn increases drastically, so is the Mn-C bond length as shown 
in Fig. Ic. Moreover, the graphene layer around the SV region becomes strongly rippled 
and two of the three carbon atoms around the vacancy site form a pair (Fig. le and 2d). 
The C dangling bond states now becomes Db, D^, and Dab (Fig. 2d). Because of the 
elongated structure, the hybridization between the Dn and Mn d orbitals becomes much 
reduced in the rippled structure so that the Mn atom maintains its d^'^d^^ configuration 
with the minimal d-level broadening for both spin-up and spin-down channels. Therefore, 
the spin transition of Mn-SV system indicates that the rippled graphene structures (induced 
by compression) play an important role in the spin transition from low-spin to high-spin 
state. It is quite encouraging to see that the rippling of graphene layers already can be 



controlled by substrates or thermal expansion in experiments 
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The interaction of localized Mn moments is evaluated by calculating the energy differ- 



TABLE I: The energy differences between ferromagnetic (FM) state and antiferroniagnetic (AFM) 
state {Eafm — ^fm) for Mn atoms on pristine graphene as a functional of strain rj. The Mn-Mn 
distance is fixed at ~ 9 A. The corresponding Curie temperature Tc are also shown. 
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ence (-Eafm — E^m) between ferromagnetic (FM) state and antiferroniagnetic (AFM) state. 
Taking Mn-PG system as an example, Table I shows the calculated -Eafm ~-Efm per Mn-Mn 
pair as a functional of r]. The Mn-Mn distance is set to ~ 9 A, which is a reasonable doping 
concentration in practice. The calculated -Eafm — -Efm can be related to the Curie temper- 
ature Tc through an effective Heisenberg model with an interaction J = (-Eafm — -Efm)/2. 
The positive energies in Table I indicate FM spin alignment at low temperature while nega- 
tive values are AFM ones. Without strain, Mn atoms on graphene prefer the AFM coupling 
in agreement with previous results^]. This is consistent with the fact that the spin-up 
d orbital is fully occupied and spin-down d orbitals are fully empty. The spin hopping is 
allowed in the AFM arrangement but not allowed in the FM configuration, resulting in a 



lower energy of AFM state [23|]. Interestingly, strain dramatically changes the spin coupling 
between Mn atoms on graphene when the Mn spin-down d orbitals become partially occu- 
pied and spin hopping is allowed for FM arrangement. For t] > 0.05, the magnetic ordering 
converts from AFM to FM, consistent with the expectation of the band coupling model 



description of the magnetic interactions 23|] . This result strongly demonstrates that the spin 
exchange strength and magnetic ordering of a TM-graphene system could be switched by 
simply using a strain. The Curie temperature, Tc, can be roughly estimated by using the 
formula Tc = 0.44J5'(5' + 1)/A;_b which correspondings to a quantum Monte-Carlo (QMC) 
study of a quasi two-dimensional Heisenberg ferromagnet on a square lattice [2J]. The results 
are shown in Table I. For example, a value of Tc = 151 K is obtained when rj = 0.05. 

Not only for the Mn atom on graphene, the strain control of the magnetism generally 
exists for almost all the TM atoms from Sc to Zn on graphene. The range of magnetic 
moments of TM-graphene systems under realistic strain -0.05 < r^ < 0.10 are shown in Fig. 3. 
The transition of magnetic states for TM-PG systems under strain could be either continuous 
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FIG. 3: The range of magnetic moments of TM atoms adsorption on graphene (a) without vacancy, 
(b) with SV, (c) with DV under -0.05 < r] < 0.10. The intrinsic magnetic moments of TM-graphene 
without strain are shown as green hne, while the maximum and minimum magnetic moments of 
TM-graphene caused by compressive (or tensible) strain (in unit of %) are shown by blue (or red) 
line. 

(e.g., V, Ti, and Cr) or abrupt (e.g., Sc, Mn, Fe, and Cu), as shown in Fig. 3a, which is 
mainly due to different strain dependent ligand splitting and hybridization. Interestingly, 
strain can also control the adsorption sites of TM atom, e.g., the adsorption sites of Cu (Sc 
and Fe) on graphene converts from C-C bridge (hollow) site to hollow (C-C bridge) site at 
7] > 0.08 (r/ < -0.04). 

The spin transition is even more remarkable in the TM-SV cases, as shown in Fig. 3b. 
For example, a compressive strain of -0.04, -0.04, and -0.035 could transfer Ti-, V- and 
Cr-SVs from low-spin to high-spin states, respectively, similar to Mn-SV. A tensile strain of 
0.025, 0.06, and 0.07 could transfer Cu-, Sc-, and Ni-SV from low-spin to high-spin states, 
respectively. It should be noticed that at least 66% (50%) of total magnetic moment of 
Cu-SV (Ni-SV) system origins from the nearest C atoms around vacancy. We have also 
studied the cases of TM-DV, since the stability of TM atoms adsorbed on DV is comparable 
to that of SV[7]. The spin transition in TM-DV system is similar to that of TM-SV but 
with a smaller critical rj for spin transition, as shown in Fig. 3c. Especially, a compressive 
strain of -0.01, -0.01, and -0.015 could convert the spin-states of Ti-, Cr-, and Mn-DV from 
low-spin to high-spin, respectively. Surprisingly, a very small tension of 0.005 could transfer 
V-DV from high-spin to low-spin state. 

In conclusion, a concept of controlling the magnetism of TM-graphene by simply using 
a strain is reported. Our calculations demonstrate that the strain-induced change of the 
local atomic structure of graphene around a TM atom makes a dramatic effect on both the 



8 



spin state of the TM atom and the exchange couphng between neighboring TM magnetic 
moments. A strong spin-dependent hybridization between TM d and graphene p orbital 
states, derived from the orbital selection of the local lattice symmetry, is responsible for the 
determination of the local electronic and magnetic structure. Our results indicate that the 
strain is available to control the magnetism of nanoscale materials for spintronics. 
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